Non-Meissner electrodynamics and knotted solitons 
in two-component superconductors 



Egor Babaev 

Physics Department, University of Massachusetts, Amherst MA 01003, USA 
Department of Theoretical Physics, The Royal Institute of Technology 10691 Stockholm, Sweden 

(Dated: February 15, 2009) 

I consider electrodynamics and the problem of knotted solitons in two-component superconduc- 
tors. Possible existence of knotted solitons in multicomponent superconductors was predicted several 
years ago. However their basic properties and stability in these systems remains an outstandingly 
difficult question both for analytical and numerical treatment. Here I propose a new perturbative 
approach to treat self-consistently all the degrees of freedom in the problem. I show that there exists 
a length scale for a Hopfion texture where the electrodynamics of a two-component superconductor 
is dominated by a self-induced Faddeev term, which is a stark contrast to the Meissner electrody- 
namics of single-component systems. I also show that at certain short length scales knotted solitons 
in two-component Ginzburg-Landau model are not described by a Faddeev-Skyrme-type model and 
are unstable. However these solitons can be stable at some intermediate length scales. I argue 
that configurations with a high topological charge may be more stable in this system than low- 
topological-charge configurations. In the second part of the paper I discuss qualitatively different 
physics of the stability of knotted solitons in a more general Ginzburg-Landau model and point out 
the physically relevant terms which enhance or suppress stability of the knotted solitons. With this 
argument it is demonstrated that the generalized Ginburg-Landau model possesses stable knotted 
solitons. 



Quintessential and unusual properties of a quantum 
fluid consisting of a large number of particles can be de- 
scribed at macroscopic length scales by a simple com- 
plex scalar field 'l'(r), or in the case of a multicom- 
ponent quantum fluid by a multiplct of complex fields 
^'i(r), i — 1,2, ..,iV. The fact that this simple descrip- 
tion is possible implies that the flow of the macroscop- 
ically large number of particles comprising a quantum 
fluid is severely constrained. The constraint in question 
is the superflow quantization condition which originates 
from the single- valuedncss of the complex fields ^'i(r). In 
a physical situation where there is a superflow locally in 
space this flow should satisfy the single-valuedness con- 
dition of the condensate wave function. Thus it can be 
created by exciting vortex loops with quantized super- 
fluid velocity circulation, or magnetic flux. Consequently 
many key properties of quantum fluids depend dramat- 
ically on properties of the vortex loops. In presently 
known and well investigated superfluids and supercon- 
ductors the main characteristic which vortex loops have 
in common is the fact that their energies depend mono- 
tonically on the loop sizes. Indeed this fact is very im- 
portant for physics of how a quantum fluid restores sym- 
metry (via entropy driven proliferation of vortex loops), 
reacts to a quench (via a relaxation of a quench induced 
vortex loops), for physics of superfluid turbulence, etc. 
All these properties would be quite dramatically altered 
if some quantum fluid would allow vortex loops with non- 
monotonic energy dependence on the loop size (i.e. if the 
energy of the vortex loop would grow not only if the 
vortex loop expands but also if it shrinks below some 
characteristic size). In these cases quench- induced de- 
fects would be protected from decay by a potential bar- 



rier resulting in large-scale remnant post-quench vortic- 
ity. It may also produce hysteretic behavior in entropy- 
generated topological defects and thus change the order 
of the superconducting phase transitions etc. Essentially 
in many respects this would lead to a new type of super- 
fluid behavior, but the outstanding question is whether 
any quantum fluids can, in principle, support such de- 
fects. 

The model which is considered below applies, with 
some modifications, to a large variety multicomponent 
systems such as multicomponent electronic condensates 
(e.g. two-band superconductors [l[), physics of the pro- 
jected multicomponent quantum metallic fluid of hydro- 
gen or its isotopes under high compression similar 
situation may arise in spin-triplet superconductors 0], 
similar models were also discussed in the context of neu- 
tron star interior [1, Q . Besides that in condensed mat- 
ter physics there is a growing interest in systems where 
SU{2) Ginzburg-Landau functional appears as an effec- 
tive model 0, 01 . In that context it is indeed important 
to understand the basic properties of topological defects 
in SU{2) superconductor in order to understand fluctu- 
ations and critical behavior in these systems. 

In 3 and also in 0, 0] it was conjectured that some 
multicomponent superconductors may support defects in 
the form of loops or knots which energy would grow if 
such a vortex shrinks. However this energy scaling and 
the question of stability of these defects in these systems 
turned out to be extremely difficult, and in spite of mul- 
tiple attempts to solve the problem no conclusive results 
were found so far. In this work I use a new approach to 
show that indeed in certain cases, condensed matter sys- 
tems should allow topological defects endowed with such 
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properties. 

The question of the existence of topological solitons 
in the form of loops or knotted loops which energy is a 
nonmonotonic function of size was first raised in math- 
ematical physics several decades ago Faddeev pro- 
posed a model consisting of a three-component unit vec- 
tor n = (711,712, Ms), |n| = 1: 



F = J d\l^^{dnf + c^eMjn-d^nxd,nf^ (1) 

where ekij is the Levi-Civita symbol. This model sup- 
ports topological defects in the form of closed or knotted 
loops (called knotted solitons) characterized by a non- 
trivial Hopf invariant For such a defect in eq. H]), 
one can expect the energy coming from the first term 
(second order in derivatives integrated over the space) to 
scale as [r] while the energy coming from the last term 
(which is fourth order in derivatives) to scale as [1/r]. 
For this reason it was argued that for this soliton there 
is an energetically preferred length scale set by the coef- 
ficient c 0, 1^, m. Numerically the solutions for these 
defects were found only a decade ago by Faddeev and 
Niemi, confirming that the energy of a knotted soliton 
[lo| has a global minimum at a certain length scale. This 
was followed by a decade of highly nontrivial numeri- 
cal and mathematical studies which uncovered a number 



of extremely interesting properties of these defects [12 1 
(for movies from numerical simulations by Hietarinta and 
Salo see ref. [l3|)- 

The outstanding question is whether such topological 
defects can be found in condensed matter systems. In 
this case several principal problems with realizability of 
these solitons can be immediately identified. First of all 
the nonmontonicity of the energy originates in a compe- 
tition between second-order and fourth-order derivative 
terms in an energy functional. In quantum fluids, the 
Ginzburg-Landau or Gross-Pitaevskii energy functionals 
are effective models where a second order gradient term 
arises from microscopic considerations via a derivative ex- 
pansion. If one would try to stabilize knotted solitons by 
obtaining higher order terms in a derivative expansion, 
the stabilization length would be such that the second 
and fourth order terms would be of the same order of 
magnitude as sixth-order and higher terms. Therefore 
the derivative expansion fails and no Ginzburg-Landau 
description exists at this length scale. Thus for a re- 
alization of the model like ([T]), the fourth-order deriva- 
tive term should have a nonperturbative origin. That is, 
there should exist a regime where second- and fourth- 
order terms have similar magnitude, at the same time 
being much larger than any higher-order terms. In Q it 
was observed that the two-component Ginzburg-Landau 
(TCGL) model can be mapped onto a model contain- 
ing the terms ([T]). There, a fourth order term in deriva- 
tives originates nonperturbatively as a contribution to 
magnetic field energy density. However in contrast to 



eq. ([T]) the fourth-order term is coupled to another field. 
The role of that complicated coupling to the additional 
massive vector field was not known. This question has 
also turned out to be very difficult to address numeri- 
cally because the studies of these solitons in the TCGL 
model are very computationally demanding. The numer- 
ical work, performed so far, explored a limited range of 
parameters [l^ (see also the early work where however 
a dimensionality-reducing axially symmetric ansatz was 
used [Til l ) . These first numerical works did not find indi- 
cations of the overlap of properties of the TCGL model 
and the Faddeev model H]), even though the results of 
the model ^ were recovered by introducing a constraint 
which straightforwardly suppresses the additional mas- 
sive vector field. On the other hand these simulations 
did not rule out that there is a parameter range where 
knotted solitons are stable in the TCGL model and the 
question remained open. In this work I address this prob- 
lem. 

Let me briefiy outline the mapping @ of the TCGL 
model onto a model containing a version of eq. ([I]). In 
the simplest form, (used in Q) the TCGL energy density 
is: 



Ji=l,2 



2A)*„p + y(|*„|) + i(VxA)2 (2) 



where = |\I'„|e*'^" are complex scalar fields which are 
coupled by the gauge field A. The symmetry breaking 
potential term V can be quite general. Its role how- 
ever is quite straightforward to evaluate. Since in what 
follows we will focus on the most interesting processes 
where the magnetic energy competes against kinetic en- 
ergy of superflow, the potential term will be used in the 
simplest SU{2) form: V = u(|*ip + |«'2p - [const]^)^ 
with a large coefficient v. In conclusion I briefly com- 
ment on the cases of the effective potentials where the 
SU{2) symmetry is broken. The equation for supercur- 
rent which follows from ^ is: 

, 1 



(3) 



Lets introduce the following notations 



-2 - |VI/J2 + 



Xn = IXnle'*"", 
\Xn\ - 

C = J/(ep2), 



Then we define the vector field 

n — (cos((/)i — 02) sin 6, cos(0i — 02) sin0, cos( 



(4) 
(5) 
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for which the following identity holds 



1 



|Vxi 



1. 



-p'iVnf. (6) 



The kinetic and magentic energy density terms of the 
model ^ then can be rewritten as 



F 



-{Vnf 



1 
2^ 



efcij ( diCj + -n • diV^ x (9jn 



V (7) 



A knotted soliton in this model is defined as a texture 
of n characterized by a nontrivial Hopf invariant. In the 
simplest case it is a "vortex loop" which can be described 
as follows. Consider a cross section of this loop. If at 
r — > oo the value of n(r) is = (n^jTijjn-j), then in 
the center of this cross section we have n'^ = — = 



In what follows I will call the point 



where n'^ = — n° = (— nj, — n^, — n") the "core", even 
thought it does not have the same meaning as the core 
of an Abrikosov vortex. Further, in this simplest case 
if we follow some path around the core where n 7^ n° 
and n n'^ a two component vector 1 (defined as a pro- 
jection of n to a plane perpendicular to ri°) winds N 
times along that path. Besides that 1 winds M times 
along any closed paths in toroidal direction, i.e. along the 
core (that means that this closed vortex is "prepared" by 
twisting a skyrmion-like fluxtube M times before gluing 
its ends). Because of these windings, for any value of n 
one can identify one or several closed helices in the phys- 
ical space where n = n (which is called a preimage of 
n). The winding N specifies the number of these helices 
while the winding M specifies how many steps these he- 
lices have. In general ease the vortex can have a form of a 
closed knot. Then the preimages of n are closed knotted 
helices (see Fig. 1). 
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FIG. 1: (Color online) A schematic picture of a toroidal knotted soliton (shown on the right). Any cross section of the flux tube 
(rectangular area, schematically shown in yellow) produces a skyrmionic texture of the three-component vector n. The texture 
of ii, close to the vortex center in the cross section, is shown on the left. The fluxtube is twisted: the red helix schematically 
shows the preimage of the north pole: i.e. where in physical space n assumes the position corresponding to the "north pole" 
on (i.e.n = (0,0,1) ) 



Consider a toroidal knot with diameter much larger than 
the penetration length A = 1/ep. There is quantized 
magnetic flux carried by this vortex. This follows from 
the following argument. Consider a vortex loops of size 
L 3> A. Consider a path tr which is much smaller than L 
such that it goes once around the core at a distance much 



larger than A . Since by definition such a path is outside 
the texture where n has windings or varies appreciably, 
we have everywhere along that path n w n° (on the left 
side of the Fig. 1 it would correspond to a contour which 
is far enough from the core so that along that contour, the 
vector n does not deviate from the value on equator in 
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the order parameter space 5^ and pointing to the right) . 
From there is follows that along this path V/3 = 
and §^'^{4>i — 02) = 0. Further if a cross section of 
the flux tube is characterized by the winding number 
N , it should contain N preimages of the north pole of 
(defined as a point in the cross section where = {] 
and thus l^"!] = 0) as well as N preimages of the south 
pole of (defined as a point in the cross section where 
6 ^ TT and thus |4'2| = 0). The only proper mapping 
from Ginzhurg- Landau variables ^i, A. to the variables 
p, n, C should be supplemented by imposing the condition 
that singlevaluedness of is preserved. This condition 
is: along any path fi in the physical space which encircles 
Ni preimages of the north pole of S'^ and N2 preimages of 
the south pole of S"^ the conditions should hold V0i = 
2ttNi and §^'V(j)2 = '2ttN2- Note that to be consistent 
with a Hopf map each of the phases should have a single 
2tt winding per zero Nn of the order parameter for 
all boundary conditions, except the case where n assumes 
positions corresponding to the north or south pole on S*^ 
at infinity and in the core. Thus using the additional 
singlevaluedness conditions we have for the path a: 



) = 0; i V(0i + 02) = 47r7V (8) 



This determines the magnetic flux enclosed by the flux- 
tube, as follows from the equations of motion Q: 



<i> = J) dsA = 
1 

_ 1 

~ e2(|*i|2 + |«-2|2) 



ds 



ds 



(9) 



where $0 = ^ is the magnetic flux quantum. 

Note that the flux quantization of a topologically non- 
trivial texture of n appears only if the model ([7| is supple- 
mented with the additional conditions of singlevaluedness 
of the original phases 0„ . Without this extra conditions 
the model ^ allows topologically nontrivial textures of 
n which do not carry magnetic flux. It is thus an impor- 
tant circumstance that the model (l7| has configurations 
which do not correspond to physical configurations in ([2]) , 
therefore the proper mapping from ([2]) to ([7]) should in- 
volve the constrain associated with the implementation 
of the singlevaluedness conditions on "ifn- Note also the 
texture carries one flux quantum per couple of spatially 
separated zeros of VPi and ^'2, therefore a knotted soli- 
ton can be interpreted as a special bound state of twisted 
fractional flux vortices [H . 

The described above topological objects are always 
well defined when p 7^ at any point in that texture. 



The more complicated question is how the energy of these 
defects scales as a function of their sizes. I.e. if these de- 
fects can minimize their energy by shrinking or if they 
will be protected from shrinking by an energy barrier. 
Lets return to the eq. ([7]). The last term is the most 
interesting here. It represents the magnetic field energy 
density which has the contributions from the massive vec- 
tor field C and also a Faddeev term [ekijH ■ diii x djU]'^ . 
It suggests that, if there are conditions where the contri- 
bution from V X C is negligible, then the magnetic field 
energy density in the TCGL model would scale as a Fad- 
deev term. However, at least in the limit e — > 0, when 
B — > one finds that ekijdiCj = —jEkij^ • diii x ^j-n 
and the model is reduced to a Gross-Pitaevskii equation 
for two decoupled fields without a self-generated Faddeev 
term. Moreover, similar behavior has been seen in the re- 
cent numerical studies of the eq. ([2]) for a range of 
parameters. This raises the question whether the elec- 
trodynamics of two-component superconductor actually 
possess a self-generated Faddeev term or if it is compen- 
sated by the field C making scaling properties of vortex 
loops monotonic. Below I develop analytic treatment for 
the role of the field C to answer this question. 

Let me first observe that if there exists a regime where 
the system tends to the model ^ (i.e. where the field 
C plays a role of small correction), then a self-consistent 
perturbative scheme for treating the field C is possible. 
The vector n indeed depends only on gauge invariant 
quantities such as the phase difference between the con- 
densates and the relative densities. Therefore one can 
always define a texture of n, and the accompanying con- 
figuration of the field C will be determined by a configu- 
ration of the vector potential corresponding to the energy 
minimum of the functional ([2|) for a given texture of n. 
Consider now a texture of n corresponding to a knot- 
ted soliton. Further consider shrinking of that texture. 
In this process the contribution from efcy-n • diH x dju 
in the energy functional will grow. The field C should 
assume then an optimal configuration from the point of 
view of the energetic balance between (i) the best pos- 
sible compensation of the efey-n - diii x djii contribution 
in the last term in ^ , and (ii) the accompanying energy 
penalty in the second term in ([7|). Observe now that 
if there exists a parameter regime and a characteristic 
size of a knotted soliton where the self-generated mag- 
netic field consists primarily of the Faddeev term contri- 
bution: Bk ~ j^ekijA ■ diii X djfi, then a contribution 
from the field C can be estimated perturbatively. Since 
it is not a priory known if such a regime exists, the per- 
turbative scheme should be self-consistent. That is, the 
breakdown of the self-consistency criterion should signal 
the violation of the assumption. Such a self-consistent 
perturbative estimate of the contribution from the field 
C can be made by using the condition that the Maxwell 
equation J = V x B should be satisfied. First observe 
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that from J = V x B it follows that: 



C = A^<^VxVxC + -Vx ekijH ■ ^^n x djU } , (10) 



and 



VxC = A^|VxVxVxC + 

i V X V x [ekijH ■ diH X djii] | 

where A = 1/ep. By substituting this into the expres- 
sion for and repeating the procedure iteratively, an 
expansion in powers of A can be generated: 



5! 

2 



^^|e/cijn • diH X djn + 

A^V X V X [efcy-n • x 3^-11] + ...}" (11) 



Applying a similar expansion for the term {p'^/2)C'^ with 
the help of eq. ()10p it follows that when A is much smaller 
than the characteristic size L of the texture of n the 
dominant terms in the model are: 

Fl^»'^ « ^(Vn)2 + ^ (,,,> . m X d,nf (12) 

On the other hand for L « A the above self-consistent 
exclusion of the field C fails, and in the limit L/X ^ 
one approaches the e = scaling where V x C « — in • 
diH X djii, as can be seen from the equation 



C f 1 

^ = <VxVxC + -V X ekijfi ■ d^ii x djU 

^0 (X^oo) 



(13) 



From here a conclusion follows that a texture character- 
ized by a nontrivial Hopf in the TCGL model with size 
much smaller than A can minimize its energy by shrink- 
ing. It means that for a given Hopf charge a knot soliton 
in ([2) can represent at most a local energy minimum, 
while the knotted solitons in the original Faddeev model 
always correspond to a global minimum (see also remark 
0). 

However a knotted soliton texture which is much larger 
than A should have energy scaling similar to that in the 
model ([!]), i.e. receiving a contribution {ekijii-dinxdjn)^ 
in the energy density. Considering the case of the lowest 
Hopf charge Q, the energy functional has stable soli- 
tons with characteristic size given by the ratio between 
the second- and fourth-order derivatives terms which is 
Lo = l/(2ep) = ^A, note however that the self-consistent 
procedure of excluding the field C, which was used to ob- 
tain breaks down at the scales of order of A, (where 
the energy of a knot texture is affected significantly by 



the coupling to C) . The lack of small parameters in this 
regime makes it difficult to estimate energy scaling. How- 
ever it is possible that there is a finite potential barrier 
at this length scale which can prevent a texture from a 
shrinkage. The observation of the instability of knotted 
solitons in the TCGL model in the recent numerical stud- 
ies [l3| may originate from the simulations being in the 
parameter range of the monotonic scaling regime. This 
will be the case if for example the initial texture is too 
small compared to A, or it can also result from A being too 
small compared to the numerical grid spacing making the 
stabilization length being too small to resolve on a nu- 
merical grid etc. Indeed finding a finite potential barrier 
in numerical simulations is a much harder task than the 
identification of the infinite energy barrier in the original 
Faddeev model (where a texture scales to the global min- 
imum no matter what are the initial conditions). Finding 
a small potential barrier with procedures like in [iJl re- 
quires fine-tuning of the initial texture to be close enough 
to the one corresponding to the local minimum. 



Importantly, the above considerations are restricted to 
solitons with lowest Hopf charges Q. However, one of 
the most remarkable properties of the model ^ is the 
existence of a Vakulenko-Kapitanskii bound, which states 
that energy of a solit_on depends on the Hopf charge Q as 
E > const ■ c|(3|3/4 



17| . Because of the 3/4 power, soli- 



tons with high Hopf charge are stable against decay into 
several solitons with lower Hopf charge. Indeed high-Q 
solitons have size which depends not only on the length 
scale given by the ratio of the coefficients in front of the 
second- and fourth-order terms in ([7]) but also on Hopf 
charge. The textures should in general be larger for larger 
Q, (though there is no simple scaling because at high-Q 
solitons develop very complicated forms [l^V It indi- 
cates a possibility that in case of a high Hopf charge, 
knotted solitions in the TCGL model may have a local 
energy minimum at the length scales larger than A, and 
thus be reasonably well described by the effective model 



The complicated nature of the energy scaling of the 
knotted soliton texture in the simplest GL model ([2]) is 
connected with the fact that it has only one length scale 
A. This scale sets the ratio between second-order and 
forth-order derivative terms in n and at the same time 
it sets the inverse mass for the vector field C. In general 
physical systems the TCGL model includes other terms 
consistent with the symmetry. One very generic term 
which is second order in derivatives is the intercomponent 
current-current interaction (the Andreev-Bashkin terms) 
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Ml. TCGL with these terms is: 



model ((21) • This follows from the following identity: 



F 



E - 

n=1.2 



(V + ieA)^'„|2 + -i/(|$^ 



-(V X 



+ a*i*2(V + ^eA)*! ■ (V + ieA)'I'2 

+ a*i*2(V - ieA)*i • (V - ieA)*2 

- a^i^-aCV - ieA)^'i • (V + ieA)1'2 

- a*;;;*2(V + ieA)*i • (V - ieA)*2 



(14) 



In physical systems the coefficient a can vary in a wide 
range and be either negative or positive. Let us con- 
sider electrodynamics and knotted solitons stability in 
the model p4|) . A new separation of variables can be 
introduced by defining the following fields: 



Xi 



u 



Xi 



X2 



sin-e^^S X2=cos-e»*% 



(Xl,X2)i? 



|*i|2(l-4a|*2|2) 



(1 - 4a|*2P)|*i|2 + (1 - 4a|*i|2)|*2P 



|*2P(l-4a|*i|2) 



' (1 - 4a|^2P)|*i|2 + (1 - 4a|*i|2)|*2P ' 
^{xiVxt - xtVxi + X2VX2 - X2VX2} + eA 
(1 - 4a|«'2ni«'iP + (1 - 4a|*i|2)|^'2p; (15) 



Lets, as in the previous example, consider the regime 
where one can neglect density fluctuations except the rel- 
ative density fluctuations described by 6 (that is, we will 
be working with the 0(3) field u coupled to a massive 
vector field C).Thcn the model can be rewritten as: 



1 
2? 



^-(VU)2 + ^C2 

V X C + -EkijU ■ diU X djU 



-4 1X2 Vxt + xtVX2|' + 1X2 Vxi - X1VX2I' . 

(l-4a|*i|2)(l-4a|*2p) ^ ' 



From here it follows that the model p4)) can be repre- 
sented as the model ([7]) with a renormalized characteris- 
tic length scale A = l/(ep) and the additional term (the 
last term in (|16|l ). This term plays a crucial role, namely 
it breaks the single-parameter character of the TCGL 



Y(Vu)^ + 

~4lX2VxI +x|Vx2|^ + 1X2 Vxi -XiVx2|^ 
ap — 



(1 -4a|^'i|2)(l -4a|*2| 



(V6l)2 + sin^ 6'(V(0i -(/.2))^ 

"2(1 -4a|^'i|2)(l -4a|^'2|2) ^ 
i±^(V0y +sin2 0(V(0i - 02))= 



(17) 



For a < the second term diminishes the energy coming 
from the term ^(Vu)^. For this reason the second-order 
derivative terms are balanced by the fourth order term 
~ {ckijU ■ diU X dju)^ at a larger texture size (relative to 
the scale A associated with the field C). Thus the model 
(fT4)) has a tunable disparity of the characteristic lengths 
for u and C. Therefore for a large enough negative a a 
knotted soliton in the model (|14p should be stabilized at 
the length scale much larger than A and thus approach 
the properties of the knotted solitons in the model ([1]). 

In conclusion, normally the essence of electrodynamics 
of a superconductor is understood as the fact that the 
phase field and vector potential combine to produce a 
massive vector field which then describes the electrody- 
namics of the system. In this work the electrodynam- 
ics of a two-component superconductor is considered and 
it is shown that a two-component superconductor not 
only has distinct electrodynamics manifested in the gen- 
eration of a Faddeev term along with a massive vector 
field, but also that there are regimes where the electro- 
dynamics is dominated by this term. There is a crossover 
to smaller length scales where effectively the electrody- 
namics does not feature the contribution in the form of 
the Faddeev term. Therefore the knotted solitons in the 
TCGL model may be a local minimum in the energy for 
a given Hopf charge (in contrast to a global minimum 
in the case of the model ([T|)). In the second part of the 
paper I showed that a more generic TCGL model with 
physically relevant mixed gradient terms possesses two 
characteristic length scales. This makes the potential 
barrier for knotted solitons tunable. Such vortex loops, 
which have a potential barrier against shrinkage, should 
lead to entirely different quench reaction, supcrfluid tur- 
bulence and physics of thermal fluctuations. The results 
may be relevant for a variety of physical systems where 
two-component Ginzburg-Landau model is realized rang- 
ing from electronic multicomponent superconductors to 
the projected mixtures of protonic and electronic con- 
densates in liquid metallic hydrogen Q • If in the models 
like the knotted solitons being a local minimum of 
the energy functional, have energy proportional to Q^^^ 
like their kin in the Faddeev model, then these defects, if 
e.g. induced by fluctuations will have tendency to pilcup. 
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Such a behavior may be relevant for understanding the 
recent observations in numerical simulations of a discon- 
tinuous phase transitions in S'C/(2)-superconductor 
In connection with applicability to physical systems, it 
should also be noted that in principle it is not necessary 
to have the exact SU{2) symmetry for the realization of 
the physics discussed above. Because a knotted soliton 
is a closed loop, it does not produce any phase windings 
at infinity, therefore an effective potential which breaks 
SU{2) symmetry to U{1) x U{1) (or even softly breaks 
it to U(l)) introduces only a finite energy penalty for a 
knotted soliton and docs not necessarily destroys it. 

I thank Juha Jaykka for the plot of the cross section 
of knot soliton. 
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